Fake plunges are very eccentric real EMRIs in disguise 
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Abstract. The capture of a compact object in a galactic nucleus by a massive black hole (MBH) is the best 
way to map space and time around it. Compact objects such as stellar black holes on a capture orbit with a 
very high eccentricity have been wrongly assumed to be lost for the system after an intense burst of radiation, 
which has been described as a "direct plunge". We prove that these very eccentric capture orbits spend actually 
a similar number of cycles in a LISA-like detector as those with lower eccentricities if the central MBH is 
spinning. Although the rates are higher for high-eccentricity EMRIs, the spin also enhances the rates of lower- 
eccentricity EMRIs. This last kind have received more attention because of the fact that high-eccentricity EMRIs 
were thought to be direct plunges and thus negligible. On the other hand, recent work on stellar dynamics has 
demonstrated that there seems to be a complot in phase space acting on these lower-eccentricity captures, since 
their rates decrease significantly by the presence of a blockade in the rate at which orbital angular momenta 
change takes place. This so-called "Schwarzschild barrier" is a result of the impact of relativistic precession on 
to the stellar potential torques, and thus it affects the enhancement on lower-eccentricity EMRIs that one would 
expect from resonant relaxation. We confirm and quantify the existence of this barrier using a statitical sample 
of 2,500 direct-summation A'-body simulations using both a post-Newtonian but also, and for the first time, 
a geodesic approximation for the relativistic orbits. The existence of the barrier prevents "traditional EMRIs" 
from approaching the central MBH, but if the central MBH is spinning the rate will be anyway dominated by 
highly-eccentric extreme-mass ratio inspirals, which insolently ignore the presence of the barrier, because they 
are driven by two-body relaxation. 



1 Introduction 

One of the most exciting results of modern astronomy is 
the discovery, mostly through high-resolution observations 
of the kinematics of stars and gas, that most, if not all, 
nearby bright galaxies harbor a dark, massive, compact ob- 
ject at their centers. 1 1 , 2 1. The most spectacular case is our 
own galaxy, the MiUcy Way. By tracking and interpreting 
the stellar dynamics at the centre of our galaxy, we have 
the most well-established evidence for the existence of a 
massive black hole (MBH). The close examination of the 
Keplerian orbits of the so-called "S-stars" (also called SO- 
stars, where the letter S stands simply for source) has re- 
vealed the nature of the central dark object located at the 
Galactic Center. By following one of them, S2 (S02), the 
mass of SgrA* was estimated to be about 3.7 x 10^ Mq 
within a volume with radius no larger than 6.25 light-hours 
13] H). More recent data based on 16 years of observations 
set the mass of the central MBH to ~ 4x10'' Mq (|5]|6l|2l[l. 
Observations of other galaxies indicate that the masses of 
MBH can reach a few billion solar masses. The existence 
of such a MBH population in the present-day universe is 
strongly supported by Soltan's argument that the average 
mass density of these MBHs agrees with expectations from 
integrated luminosity of quasars (|9] [TOl . 
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To interact closely with the central MBH, stars have 
to find themselves on "loss-cone" orbits, which are orbits 
elongated enough to have a very close-in pericenter ifTTl 
[121 [13 1 . The rate of tidal disruptions can be established 
(semi-)analytically if the phase space distribution of stars 
around the MBH is known IIT4l [TSl [T6l for estimates in 
models of observed nearby nuclei. To account for the com- 
plex influence of mass segregation, collisions and the evo- 
lution of the nucleus over billions of years, detailed numer- 
ical simulations are required, however Iil7i il8i ,19. i20i ,211 

Many correlations linking the MBH's mass and overall 
properties of its host spheroid (bulge or elliptical galaxy) 
have been discovered. The tightest are with the spheroid 
mass f26\, its velocity dispersion (M-cr relation, |27|) and 
degree of concentration [28 1. Consequently, understanding 
the origin and evolution of these MBHs necessitates their 
study in the context of their surrounding stellar systems. 

These observations are difficult, specially for low-mass 
MBH (ranging between 10^ and 10^ Mq) and more so for 
the even less-massive intermediate MBHs (IMBHs). Nowa- 
days using adaptive optics we could optimistically hope to 
get a handful of measurements of stellar velocities of such 
targets about ~ 5 kpc away in some ten years. Neverthe- 
less, we need a bright reference star to guide us; the re- 
quirement of an astrometric reference system is crucial. 

Therefore, in order to follow the same scheme of detec- 
tion of inner stellar kinematics and of the masses of MBHs 
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in distant galaxies, we need the Very Large Telescope in- 
terferometer and one of the next-generation instruments, 
either the VSI or GRAVITY |29 30|. Only then can we 
improve our astrometric accuracy by the necessary factor 
of ~ 10. This could allow us to follow the motion of stars 
orbiting a central MBH in a galaxy with the resolution we 
need to extract the mass of the black hole from the kine- 
matics of the stars. 

This is what makes an experiment such as eLISA 11311 
|32l so appealing from the viewpoint of an astrophysicist, 
since it can be envisaged as a magnifying glass that will 
look much deeper and with more detail in the areas of in- 
terest. However, this is only true for compact stars. Whilst 
main-sequence stars are tidally disrupted when approach- 
ing the central MBH, compact objects (stellar black holes, 
neutron stars, and white dwarfs) slowly spiral into the MBH 
and are swallowed after some ~ 10^ orbits in the eLISA 
band. At the closest approach to the MBH, the system emits 
a burst of GWs which contains information about the space- 
time geometry, the masses of the system, and the spins of 
the MBH. We can regard each such burst as a snapshot of 
the system. This is what makes EMRIs so appealing: a set 
of ~ 10'' bursts of GWs radiated by one system will tell us 
with the utmost accuracy about the system itself, it will test 
general relativity ll33l . it will tell us about the distribution 
of dark objects in galactic nuclei and globular clusters and, 
thus, we will have a new understanding of the physics of 
the process. Besides, new phenomena, unknown and unan- 
ticipated, are likely to be discovered. 

On the other hand, this process is also directly linked 
to the growth of MBHs: Although there is an emerging 
consensus about the growth of large-mass MBHs thanks to 
Soltan's argument, MBHs with masses up to 10^ Mq, such 
as our own MBH in the Galactic Centre (with a mass of 
~ 4 X 10^ Mq), are enigmatic. There are many different ar- 
guments to explain their masses: accretion of multiple stars 
from arbitrary directions 1.34,1 , mergers of compact objects 
such as stellar-mass black holes and neutron stars 1351 or 
IMBHs falling on to the MBH f36|, or by more peculiar 
means such as accretion of dark matter |37| or collapse of 
supermassive stars l38l EH |40l [41] |42l . Low -mass MBHs 
and, thus, the early growth of all MBHs, remain a conun- 
drum. 

eLISA is in this regard also attractive, for it will scruti- 
nize exactly the range of masses fundamental to the under- 
standing of the origin and growth of supermassive black 
holes. By extracting the information encoded in the GWs 
of this scenario, we can determine the mass of the central 
MBH with a ridiculous relative precision of ~ 10 Ad- 
ditionally, the mass of the compact object which falls into 
the MBH and the eccentricity of the orbit will be recovered 
from the gravitational radiation with a fractional accuracy 
of also ~ 10 All this means that eLISA will not be "just" 
the ultimate test of general relativity, but an exquisite probe 
of the spins and range of masses of interest for theoretical 
and observational astrophysics and cosmology. 

2 Spinning lABHs 

If the central MBH has a mass larger than 10^ Mq, then the 
signal of an inspiraling stellar black hole, even in its last 
stable orbit (LSO) will have a frequency too low for detec- 
tion. On the other hand, if it is less massive than IO'^Mq, 



the signal will also be quite weak unless the source is very 
close. This is why one usually assumes that the mass range 
of MBHs of interest in the search of EMRIs for eLISA is 
between [10^, 10"*] Mq. Nonetheless, if the MBH is rotat- 
ing rapidly, then even if it has a mass larger than 10^ Mq, 
the LSO will be closer to the MBH and thus, even at a 
higher frequency, the system should be detectable. This 
would push the total mass to a few ~ 10^ Mq. 

For a binary of an MBH and a stellar black hole to be 
in the eLISA band, it has to have a frequency of between 
roughly 10 and 1 Hz. The emission of GWs is more effi- 
cient as they approach the LSO, so that eLISA will detect 
the sources when they are close to the LSO line. The to- 
tal mass required to observe systems with frequencies be- 
tween 10 and 1 Hz and is of 10"* - 10^ Mq. For masses 
larger than 10^ Mq the frequencies close to the LSO will 
be too low, so that their detection will be very difficult. On 
the other hand, for a total mass of less than lO^* Mq we 
could in principal detect them at an early stage, but then 
the amplitude of the GWs would be rather low. 

As the star spirals down towards the MBH, it has many 
opportunities to be deflected back by two-body encounters 
onto a "safer orbit" |43, 44|, hence even the definition of 
a loss cone is not straightforward. Once again, the prob- 
lem is compounded by the effects of mass segregation and 
resonant relaxation, to mention two main complications. 
As a result, considerable uncertainties are attached to the 
(semi-)analytical predictions of capture rates and orbital 
parameters of EMRIs. 

Naively one could assume that the inspiral time is dom- 
inated by GW emission and that if this is shorter than a 
Hubble time, the compact object will become an EMRI. 
This is wrong, because one has to take into account the 
relaxation of the stellar system. Whilst it certainly can in- 
crease the eccentricity of the compact object, it can also 
perturb the orbit and circularize it, so that the required time 
to inspiral in, f^^, becomes larger than a Hubble time. The 
condition for the small compact object to be an EMRI is 
that it is on an orbit for which f^,^ <sc (1 - e)fj Ii44i . 
with fj. the local relaxation time. When the binary has a 
semi-major axis for which the condition is not fulfilled, 
the small compact object will have to be already on a so- 
called "plunging orbit", with e > Cpiunge = 1 - ^/Jg^j^^/a, 
where R^^,^^ is the Schwarzschild radius of the MBH, i.e. 
^Schw - 2GMJc^, with M, being the MBH mass. It has 
been claimed a number of times by different authors that 
this would result in a too short burst of gravitational radi- 
ation which could only be detected if it was originated in 
our own Galactic Center |45 1 because one needs a coher- 
ent integration of some few thousands repeated passages 
through the periapsis in the eLISA bandwidth. 

Therefore, such "plunging" objects would then be lost 
for the GW signal, since they would be plunging "directly" 
through the horizon of the MBH and only a final burst of 
GWs would be emitted, and such burst would be (i) very 
difficult to recover, since the very short signal would be 
buried in a sea of instrumental and confusion noise and (ii) 
the information contained in the signal would be practi- 
cally nil. There has been some work on the detectability of 
such bursts Il45ll46ll47ll . but they would only be detectable 
in our galaxy or in the close neighborhood, but the event 
rates are rather low, even in the most optimistic scenarios. 
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In figures [T] and |2] we show plots of the location of the 
LSO in the plane a (pc) - (1-e), including the Schwarzschild 
separatrix between stable and unstable orbits, p-6-2e - 0, 
for both prograde and retrograde orbits and for different 
values of the inclination t. Each plot corresponds to a dif- 
ferent value of the spin, showing how increasing the spin 
makes a difference in shifting the location of the separa- 
trix between stable and unstable orbits, pushing prograde 
orbits near GM, /c^ while retrograde orbits are pushed out 
towards 9GM,/c^. The procedure we have used to build 
these plots is a standard one. Briefly, given a value of the 
dimensionless spin parameter s = a,c^ /(GM,) and a value 
of the eccentricity e and inclination angle t, and we have 
followed the algorithm, definitions and notation introduced 

inia. 

In Il48ll it was estimated the number of cycles that cer- 
tain EMRI orbital configurations that were thought to be 
plunging orbits (or orbits with no sufficient cycles) in the 
case of non-spinning MBHs can spend in a frequency regime 
of / e [lO"'*, 1] Hz during their last year(s) of inspiral be- 
fore plunging into the MBH. This is important to assess 
how many of these EMRIs will have sufficient Signal-to- 
Noise Ratio (SNR) to be detectable. It was found that (pro- 
grade) EMRIs that are in a "plunge" orbit actually spend a 
significant number of cycles, more than sufficient to be de- 
tectable with good SNR. The number of cycles has been 
associated with A^^ (the number of times that the azimuthal 
angle (p advances 2n) which is usual for binary systems. 
However, as we have discussed above the structure of the 
waveforms from EMRIs is quite rich since they contain 
harmonics of three different frequencies. Therefore the wave- 
forms have cycles associated with the three frequencies 
(fr'fg'fif) which makes them quite complex and in prin- 
ciple this is good for detectability (assuming we have the 
correct waveform templates). Moreover, these cycles hap- 
pen just before plunge and take place in the strong field re- 
gion very near the MBH horizon. Then, these cycles should 
contribute more to the SNR than cycles taking place farther 
away from the MBH horizon. 

We also estimate the impact on the event rates. Since 
"direct plunges" are actually diguesed EMRIs, although 
with a higher eccentricity. We prove that 
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In the expression 'W is a function that depends on t, the 
inclination of the EMRI and s, its spirQ We also have as- 
sumed that the SBHs distribute around the central MBH 
following a power-law cusp of exponent y, i.e. that the 
density profile follows p cc r^y within the region where 
the gravity of the MBH dominates the gravity of the stars, 
with y ranging between 1.75 and 2 for the heavy stellar 
components 1 24l l25l [49l l50l |5T1 l52l l53l B4l and see ll55l 
for an interesting first idea of this conceprl 



' For the derivation and some examples of values for "W, we 
refer the reader to the work of 1481 . 

^ The authors obtained a similar solution for how electrons dis- 
tribute around a positively charged Coulomb centre. 



For instance, for a spin of s = 0.999 and an inclina- 
tion of L - 0.4 rad, we estimate that ~ 0.26 and, thus, 
^EMRi ~ boost the event rates by a factor of 30 

in comparison to a non-rotating MBH. 



3 Analysing a blockade In angular 
momentum for low-eccentricity EI\/IRIs with 
a direct-summation A^-body sample of 
2,500 simulations 

In a gravitational potential with a high degree of symmetry, 
a test star will receive gravitational tugs from the rest of the 
field stars which are not totally arbitrary and hence do not 
add up in a random walk way, but coherently. The potential 
will prevent stellar orbits from evolving in an eiTatic way. 
In a two-body Keplerian system, a stellar-mass black hole 
will orbit around the MBH in a fixed ellipse. The stellar 
BH will not feel random gravitational tugs. It evolves co- 
herently as the result of the action of the gravitational po- 
tential. When an EMRI approaches the periapsis of its or- 
bit, we can envisage the situation as a pure two-body prob- 
lem; initially Newtonian but later GR effects must be taken 
into account as the periapsis grows smaller and smaller. 
Nonetheless, as the stellar BH goes back to the apoapsis, 
it will feel the surrounding stellar system, distributed in 
the shape of a cusp which grows in mass the further away 
we are from the periapsis. The time spent in the region in 
which we can regard this as a two-body problem is much 
shorter than the time in which the stellar BH will feel the 
rest of the stellar system. This is particularly true for the 
kind of objects of our interest, since the very high eccen- 
tricity implies a large semi-major axis. The time spent on 
periapsis is negligible as compared with the time spent on 
apoapsis. In that region, the stellar BH feels the graininess 
of the potential. The gravitational tugs from other stars 
wiU alter its orbit. The mean free path in /-space of that 
test stellar BH is very large and thus, it has a/aif random 
walk. Both the magnitude and direction of J of the stellar 
black hole are altered. When the modulus is changed but 
not the direction, we talk of "scalar" resonant relaxation, 
and correspondingly when the direction is changed but not 
the modulus, "vector" resonant relaxation. 

In particular, in the potential of a point mass, orbits are 
frozen fixed ellipses that exert a continuous torque on the 
test star. A test star does not feel random kicks from all 
directions. When we add up the individual contributions 
coming from all the rest of the stars on the test star, there 
is a residual, non-negligible torque that will influence its 
evolution. The mean free path of the star in J space is very 
large. We will refer to this phenomenon as scalar resonant 
(or coherent) relaxation, because it can change both the 
magnitude of J and the inclination of the orbital plane of 
the test star In this scenario it is possible to alter an initially 
very circular orbit and modify it in such a way that the test 
star will get very close to the MBH after the torques have 
acted, i.e., we open a new window for stars to fall into a 
capture orbit that will lead to an EMRI. 

The impact of coherent relaxation on the production 
of EMRIs is important. While the underlying physics of 
the process is very robust, it is a rather difficult task to as- 
sign values to the different parameters on which the pro- 
cess depends. A possible way of assigning these was sug- 
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Fig. 1. LSO for a MBH of mass 4 x 10^ and a SBH of mass m. = 10 for a Kerr MBH of spin s = 0.4 (left) and s = 0.7 
(right). Tiie Schwarzschild separatrix is given as a solid black line. Curves above it correspond to retrograde orbits and inclinations of 
I = 5.72, 22.91, 40.10, 57.29, 74.48 and 89.95° starting from the last value (89.95°). In the left panel we can barely see any difference 
from the different inclinations due to the low value of the spin. 
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Fig. 2. As in figure[l]but for a spin of s = 0.99 (left) and s = 0.999 (right panel). The larger the spin, the "further away" the Kerr LSO 
gets from the Schwarzschild LSO. 



gested in I56ll57l . In Figure 6 of [57], they show the rate 
of EMRIs and plunges in a system in which we take into 
account both orthodox or regular relaxation and coherent 
relaxation, normalised to what one can expect when only 
taking into account normal relaxation. The numerical sim- 
ulations of 1 57 1 show that coherent relaxation can enhance 
the EMRI rate by a factor of a few over the rates predicted 
assuming only slow stochastic two-body relaxation. 

Recently, |58| estimated with a few direct-summation 
A^-body simulations expanded with a statistical Monte- 



Carlo study that the production of "traditional EMRIs" via 
resonant relaxation is markedly decreased by the presence 
of a blockade in the rate at which orbital angular momenta 
change takes place. This so-called "Schwarzschild barrier" 
is a result of the impact of relativistic precession on to 
the stellar potential torques. Although the authors find that 
some particles can penetrate the barrier, EMRIs are signif- 
icantly supressed in this scenario. 

We have felt motivated by these results and recently 
decided to study the effect with a set of some 2,500 direct- 
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summation A^-body simulations, including post-Newtonian 
corrections and also, for the first time, the implementation 
of a solver of geodesic equations in the same code |59|. 
For this, we use a modification of the publicly available 
planet code by Sverre Aarseth, a direct summation A^- 
Body integrator 16011611 . For our study we use several dif- 
ferent methods to account for the general relativistic cor- 
rections ot the Newtonian accelerations: 

- purely Newtonian 

- post-Newtonian (PN) corrections 

- geodesic equations 

In the purely Newtonian case, the integration is done 
without modifications to the acceleration equations. In the 
PN case we add PN corrections in the following way: 

Newt periapsis shift energy loss neglected 

F ^^~F\r + c'^F2 + c-^F4+ c-^Fs +0{c-^), (3) 

IPN 2PN 2.5PN 

i.e., as described in the pioneering work of 1*621. The indi- 
vidual Fi can be found in |63 1, their equation (7.16). 

Given the high mass ratio for EMRIs, their motion around 
a MBH can also be approximated by solving the geodesic 
equations. These equations give the exact trajectory of a 
test mass particle around a Schwarzschild BH. Unlike the 
PN approximation, the geodesic equations are valid even in 
the last few during a plunge or inspiral, however only in 
the limit mi,/M, — » 0. Some orbits are expected to migrate 
towards plunge or inspiral orbits at pericenter distances of 

< 15 rg, where the errors of the PN approximation can 
already be quite significant [64J. In order to test the exis- 
tence of the Schwarzschild barrier at small distances, we 
implemented these corrections in the planet code. 

Although we confirm the blockade, we note that the 
systems the two groups study are representative but artifi- 
cial, because they integrate only 50 particles with the same 
mass distributed around the MBH. In this regard, the exis- 
tence of this quenching must be still explored in more de- 
tail, in particular with direct-summation simulations with 
realistic number of stars and mass ranges. 

In any case, in this scenario the barrier can roughly be 
fitted ma- e space as 

flSB ^CsB-(l-ey'^'. (4) 

In the last equation osb is the semi-major axis in mpc be- 
low which penetration of particles is severely supressed 
and CsB is a constant of order unity. |[58l choose Csb = 0.7. 
We depict the barrier in figure[6] As we can see, and as de- 
scribed in the work of |58 1, the barrier poses a real problem 
for stars with small semi-major axis, and below it, the evo- 
lution is dominated by two-body relaxation. 

In order to quantify the nature of the "Schwarzschild 
barrier", we first plot the normalized presence density as 
a histogram in the {a, 1 - e) plane for the Newtonian case, 
Fig.|4](left panel) and the relativistic case (right panel), and 
we give the theoretical distribution in figure [3] If we con- 
sider our specific setup, there are 3 different regions in the 
{a, \ - e) plane where different mechanisms are efficient. 
In the right-most region, where pericenters are large, RR 
plays the dominant role. The left border of this region is 
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Fig. 3. Theoretical distribution for the density of presence of a 
cusp of power-law 1.75 around a MBH, using a truncated distri- 
bution. 

roughly given by the appearance of the Schwarzschild pre- 
cession which inhibits the BHs from experiencing coherent 
torques 1591. We derive that Csb = 0.35. 



4 Only fake plunges survive: Real EMRIs 

For a fast spinning MBH, the separatrix for prograde or- 
bits is shifted to significantly lower a values, with a cor- 
responding higher value of the critical semi-major axis, 
corresponding to the point PP in figure [6] As we have ex- 
plained above, it is this effect which can lead to a signifi- 
cant increase in the EMRI rate, combined with the fact that 
the critical point for retrograde orbits (PR) is much less af- 
fected and that an isotropic orbit distribution is expected, 
thanks to relaxational processes. However this increase in 
EMRI rate would can be thwarted by vector RR if this pro- 
cess can change the orbital orientation of a SBH after it has 
crossed the "f^^ - p^^j" line and before it has completed 
its GW-driven inspiral, i.e. on a timescale shorter than t^^. 
Indeed, if the orbit becomes significantly less prograde as 
the the inspiral takes place, due to RR, the separatrix moves 
up and the SBH might suddenly find itself on a plunge or- 
bit. 

To check for this possibility, we also plot, in figure |6] 
a long-dashed line coiTesponding to the condition t^^ - 

W V' ^^'^^ ^Gw ^ ^RR V °" ^^^^ 'S^^Hs that 

cross the "fr-,,, - t " line while on the left side of the 

(jW r, pen 

"^Gw ~ ^RR v" ^^^^ keep their orbital orientation during 
their inspiral and complete it without abrupt plunge. One 
can see that, for our choice of parameters, this is the case 
for all prograde orbits. On the other hand, retrograde orbits 
can cross the "fr.,,, = t " line while RR is still effective 

GW r, pen 

enough to change their orientation during inspiral. How- 
ever, the effect of RR on retrograde orbits cannot reduce 
significantly the total EMRI rate and may even increase 
it slightly because (1) these orbits contribute less that the 
prograde ones (and more to the plunge rate) and (2) sta- 
tistically, RR is more likely to make the orbit become less 
retrograde which pushes down the separatrix. 

Moreover, as we can see in figure |6] for semi-major 
axis with values approximately a ^ 0.03 pc, the barrier lies 
well below the last separatrix, and so EMRIs originating 
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Newtonian Dynamics 



Dynamics with IPN, 2PN and 2.5PN 





Fig. 4. Integrated presence density for the Newtonian (left panel) and the relativistic case (right panel). The shaded box marks the region 
of the slice analyzed in Fig. |5] The lines indicate the position of the Schwarzschild barrier with Csb = 0.35 (blue) and the limit for 
capture onto inspiral orbits for non-resonant relaxation (green). 
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Fig. 6. Left panel: Relation between different timescales in the s = 0.7 case. We display the Schwarzschild separatrix as a solid, black 
line and the separatrices for different inclinations with different curves in light grey. The dashed, blue line shows the value of a and 1-e 
for which the vectorial resonant relaxation timescale (/j^j^ ^) is equal to the gravitational loss timescale (t^^). The dashed, dotted line 
corresponds to the values of a and 1-e for which the relaxation time at periapsis (f, ^^^j) equals the gravitational loss timescale. The 
dashed-dotted, black curve describes the "Schwarzschild barrier". Right panel: Same as the left panel but for a spin value of s = 0.999. 
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Fig. 5. Integrated presence density for a semi-major axis slice 1.5 
mpc < a < 2.0 mpc. The shaded boxes represent the position of 
the Schwarzschild barrier with Csb = 0.35 (blue) and the limit for 
capture on to inspiral orbits for non-resonant relaxation (green). 



from these area, the "plunges" we have discussed in this 
work, will not suffer this quenching in phase space. 



5 Conclusions 

The event rate of "plunges" is much larger than that of EM- 
RIs, as a number of different studies by different authors 
using different methods find. Up to now spin effects of the 
central MBH have been always ignored. Hence, the ques- 
tion arises, whether a plunge is really a plunge when the 
central MBH is spinning. This consideration has been so 
far always ignored. 

So as to estimate EMRI event rates, one needs to know 
whether the orbital configuration of the compact object is 
stable or not, because this is the kernel of the difference 
between an EMRI and a plunge. In this paper we take into 
account the fact that the spin makes the LSO to be much 
closer to the horizon in the case of prograde orbits but it 
pushes it away for retrograde orbits. Since the modifica- 
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tions introduced by the spin are not symmetric with re- 
spect to the non-spinning case, and they are more dramatic 
for prograde orbits, we prove that the inclusion of spin in- 
creases the number of EMRI events by a significant factor. 
The exact factor of this enhancement depends on the spin, 
but the effect is already quite important for spins around 
s ~ 0.7. 

We also prove that these fake plunges, "our" EMRIs, 
do spend enough cycles inside the band of eLISA to be 
detectable, i.e. they are to be envisaged as typical EMRIs. 
We note here that whilst it is true that EMRIs very near 
the new separatrix shifted by the spin effect will probably 
contribute not enough cycles to be detected, it is equally 
true for the old separatrix (Schwarzschild, without spin). 
In this sense, we find that the spin increases genericaUy 
the number of cycles inside the band for prograde EMRIs 
in such a way that EMRIs very near to the non-spin separa- 
trix, which contributed few cycles, become detectable EM- 
RIs. In summary, spin increases the area, in configuration 
space of detectable EMRIs. We predict thus that EMRIs 
will be highly dominated by prograde orbits. 

We also demonstrate that these new kind of EMRIs we 
describe here originate in a region of phase-space such that 
they will be ignorant of the Schwarschild barrier. The rea- 
son for that is that they are driven by two-body relaxation 
and not resonant relaxation. While the boost in EMRI rates 
due to resonant relaxation is affected by the Schwarzschild 
barrier, so that "standard" EMRIs run into the problem 
of having to find a way to cross it, our "plunge-EMRJs" 
are already in the right place and led by two-body relax- 
ation. The barrier affects the productions of EMRIs via 
torques, but not two-body relaxation, which is the mech- 
anism producing the "plunge-EMRIs". Moreover, because 
the supression in the rates is severe for those EMRIs with 
semi-major axis with values approximately a ^ 0.03 pc, 
we predict that the rates will be dominated by the kind of 
EMRIs we have described in this work. 
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